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Abstract. - We argue that for complete wetting at a curved substrate (wall) the wall- fluid 
surface tension is non-analytic in R^, the curvature of the wall and that the density profile 
of the fluid near the wall acquires a contribution proportional to the gas-liquid surface tension 
xi?^ 1 plus higher-order contributions which are non-analytic in R~ x . These predictions are 
confirmed by results of density functional calculations for the square-well model of a liquid ad- 
sorbed on a hard sphere and on a hard cylinder where complete wetting by gas (drying) occurs. 
The implications of our results for the solvation of big solvophobic particles are discussed. 



Understanding the adsorption of fluids at solid substrates has taken on new importance 
with recent advances in the controlled fabrication of tailored surfaces for applications in mi- 
crofluidics and other areas pQ. Much experimental [2] and theoretical effort 3 is concerned 
with wetting and associated interfacial transitions in wedge geometry and recently attention 
has turned to wetting at an apex 0]. It is becoming increasingly clear that substrate geom- 
etry can have a profound influence on the nature of fluid adsorption and, in particular, on 
wetting characteristics making these quite different from those at a planar substrate. Here 
we consider complete wetting at two substrates that have simple geometries, namely a single 
sphere of radius R s and an infinitely long cylinder of radius R c . Using an effective interfacial 
Hamiltonian approach 5 combined with exact microscopic sum-rules for the density profile 
of the fluid near a hard wall, we show that in the limit Ri — > oo, with i — s, c, the surface 
tension of the substrate (wall)-fluid interface is non-analytic in the curvature i?" 1 and that 
the density of the fluid in contact with the hard wall acquires a contribution proportional 
to 7 ff ;(oo) / Ri, where 7 g z(oo) is the surface tension of the planar interface between coexisting 
gas and liquid, as well as higher-order non-analytic terms in RJ 1 ■ Our results, which are 
confirmed fully by the results of microscopic density functional (DFT) [H] calculations, show 
that non-zero curvature leads to unexpected and subtle effects on interfacial properties when 
complete wetting occurs, even for the simplest of substrate geometries. 

In previous theoretical studies of wetting on spheres and cylinders [7181911(7] the thrust was 
on understanding how a finite radius limits the thickness of a wetting film and modifies the 
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wetting transition that can occur at a planar substrate. Little attention was paid to the effect 
of curvature on the surface tension and on the form of the density profile near the substrate 
which are the main concerns of this Letter. Our results have repercussions for the general 
theory of solvation of big solvophobic solute particles, for wetting of colloidal particles and 
fibers [TlllOj and, possibly, for the surface tension of drops and bubbles 

We consider first the general case of a bulk fluid phase a, at chemical potential fi, in contact 
with a wall w and implement a standard, coarse-grained, effective interfacial Hamiltonian 
approach 0|Ej in which the wetting film of fluid phase b is characterized only by its thickness 
I. For short-ranged (finite ranged, exponentially or faster decaying) wall-fluid and fluid-fluid 
potentials the binding potential, i.e. the excess (over bulk) grand potential per unit area, 
takes the form [5J|S]: 



£U0 = iM) + ll(Ri + ~a{T)e- m + (Pa ~ Pb)S(i I + tt ' 7ta ' ' + ' + Oil/Ikf, (1) 



with i — s, c and a s = 2, a c — 1 for the spheres and cylinders, respectively. 7^ 6 is the surface 
tension of the wall-phase b interface, 7^ a is the tension of the ba fluid-fluid interface located 
near R4 + I, a(T) > is a coefficient that we need not specify further for complete wetting, 
and £ is the true bulk correlation length of the wetting phase b. The fourth term in eq. Q 
accounts for the increase in grand potential per unit area associated with the volume of a film 
of phase b; Sfj, = [i — u co (T) is the chemical potential deviation and p a , pb are the number 
densities of bulk phases a and b at coexistence /i co (T). For Ri — 00 the binding potential 
Q reduces to that appropriate for complete wetting at a planar interface in models where 
the range of the wall-fluid potential is shorter than £ and, henceforward, we shall assume this 
is the case. Then minimizing eq. Q w.r.t. I yields the well-known logarithmic divergence 
of the equilibrium film thickness: l eq (oo) = — £ln(<5/i £(p a ~ Pb)/a{T)), as S/j, — > 0. When 
Ri ^ 00 two important modifications arise: i) 7^ b and 7 ba now depend on curvature and ii) 
the surface area of the ba interface now depends on the film thickness I. Equation JQ) assumes 
very large radii Ri and that l/Ri ^ 1. Note that the fifth term in is proportional to the 
Laplace pressure across the fluid-fluid ba interface. Since 7^ is always non-zero, away from 
the bulk critical point, this term ensures that the film thickness remains finite even at /i co (T), 
i.e. minimizing eq. £Q| yields l % e JRi) = — (,\n[ai^l a (Ri)/(Ria(T))] for Sfi — 0, where we 
have ignored terms 0(l/Ri) 2 and the / dependence in 7^ - this can be justified a posteriori. 
Several authors HIEHHIEII have argued that for large Ri the quantity ai~fl a (Ri) / Ri should 
play the same role in complete wetting at a curved wall, with <5/i ~ 0, as the effective bulk field 
(Pa — Pb)$P plays for the planar wall. We shall pursue this argument further in the present 
Letter. 

The wall-fluid surface tension is given by 



7 URi) = "ia(le q ) = lURi) + llaiRi) + « + ^) + (Pa - Pb)Sp . (2) 



For the planar interface the final term vanishes at <5/i = and the wall-fluid tension reduces to 
the sum of the planar tensions: j wa (oo) = "f w b(oo) + 7b a (oo), appropriate to complete wetting 
by a macroscopic film of phase b. When <5/j =/= 0, "f W a(oo) acquires a non-analytic contribution 
proportional to — |<5/i| In \Sp\ For finite Ri we set Sfj, = and obtain 





Ri 



\n(R l x const) + H.O.T.. (3) 
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That the surface tension contains a contribution which is non-analytic in the curvature is 
clearly a direct manifestation of complete wetting. However, the necessity for such a contribu- 
tion does not seem to have been widely recognized. An exception is Ref. [7] but in that paper 
the consequences were not discussed. In order to elucidate these we must examine the other 
terms in eq. 0- At the non-wet wb interface we do not expect non-analyticities in J^^Ri), 
provided the system exhibits short-ranged forces jTJJ- For a spherical fluid- fluid interface it is 
usually assumed m (for short-ranged forces) that 

Y ba (R s ) = 7 ba(oo)(l - 26 S T /R S + H.O.T.), (4) 

where is the Tolman length familiar in studies of liquid drops, whereas for cylindrical 
interfaces it has been argued [TQ that 

iL(Rc) = 7&aM(l + b H \nR c /R c + H.O.T.). (5) 

If the latter form were correct eq. Q would imply that the wall- fluid tension 7^, (i? c ) should 
have a contribution 7fc a (oo)(frff + £) hi(i? c )/i? c . Provided the length bn is comparable with 
the bulk correlation length £, as is expected on physical grounds, the term in bn should be 
easily identified in numerical work. 

The result J2J is valid for any complete wetting situation, with the proviso that the forces 
are short-ranged. We specialize now to the case of a hard wall exerting a purely repulsive 
potential on the fluid: Vi{r) = oo for radius r < Ri and is zero for r > Ri. It is well-known 
that the phenomenon of complete drying occurs at a planar hard wall: the interface between 
the bulk liquid (I) and the wall is wet by a macroscopic film of gas (<?) as p — > pt {T) for all 
temperatures T at which gas and liquid coexist |15H16| . At p = p£ (T) the density profile 
of the fluid is a composite of the planar wall-gas and the (free) gas-liquid interfacial profiles 
and 7,^(00) = 7 u , g (oo) + 7 ff ;(oo); phase a = I and phase b = g. Whilst drying at a planar 
hard wall has been investigated extensively using computer simulations |15| and DFT 
we are not aware of studies on curved substrates. We shall see below that curvature leads to 
some striking new results. When Ri ^ 00 thick films of gas will still develop at the wall but 
the thickness will remain finite at pJ a (T) and the wall- liquid tension 7* i(Ri) should exhibit 
the non-analyticities described above. Moreover the density profile p(r) is no longer a perfect 
composite of the two separate interfacial profiles and exhibits interesting features near the 
wall. 

A particular advantage of hard walls is that there is an exact statistical mechanical sum- 
rule 53 which relates the fluid density at contact, p l (Rf) 1 to the pressure p of the bulk 
(reservoir) fluid and the wall-fluid tension: 

In the limit Ri — 00 the contact density reduces to p/fc^T, the well-known planar contact 
theorem. Equation JBJ is valid for any one-component fluid at a hard wall. We now insert our 
result {S} for the tension of the 'dry' interface, 7^, ; (i?i), into © to obtain the contact density 
Pl iq (Rf) at fi = fif Q (T). By subtracting p l gas (Rf ), the contact density when the bulk phase is 
gas at p = p,~ a (T), we eliminate both the pressure p(p co ) and the (analytic) wall-gas tension 
obtaining for a hard spherical wall: 

kBT(p] iq (Rt) - P s gas (Rt)) = + lnR s + H.O.T. (7) 
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Fig. 1 - The density profile of a square-well fluid adsorbed at a hard wall for ksT/e = 1. Solid 
curve refers to a liquid at bulk coexistence pto(T) near a spherical wall of radius R s = 2500<7hs and 
the dotted curve to the gas at u^ a (T) at the same wall. Symbols refer to the liquid near a planar 
wall at chemical potential fi — fi co (T) = 2"/ g i(oo) /(R s (pi — p 9 )) (see text). One observes that the two 
wall-liquid profiles are indistinguishable in the region of the gas-liquid interface and close to the wall 
- see inset. Contact occurs at distance h — aHs/2. 

and for a hard cylindrical wall: 

k B T{pURt) - PUW = ^ + (g + ^?' (00) + H.O.T. (8) 

where H.O.T. refers to terms of higher order in Rf . For a planar wall the contact density 
is p(fJ- C o)/kBT in both phases so the difference vanishes identically. There are two remarkable 
features in these results. First the contact density in the liquid phase ('dry' interface) depends 
on 7 g z(oo), the tension of the planar gas- liquid interface, which can be very far from the wall 
as Ri — > oo: the first term on the r.h.s. of 10 and JHJ is simply the Laplace pressure across 
the gas-liquid interface. Second the contact density difference for the sphere contains a next 
to leading order Rj 2 \nR s non-analyticity, whereas for the cylinder the next to leading order 
term is analytic, i.e. 0(R~ 2 ), despite the fact that both surface tensions are non-analytic at 
order i?" 1 lni?^ 

In order to test the predictions of the coarse-grained theory we have adopted a fully 
microscopic DFT approach and performed calculations for a square-well fluid at hard curved 
walls. The fluid-fluid pair potential <fr(r) is infinite for r < cjhs, the hard sphere diameter, 
the width of the well is 0.5ans and the well-depth is e. We treat the hard-sphere part of 
the free-energy functional by means of Rosenfeld's fundamental measure theory JSj and the 
attractive part using a simple mean- field approximation 6,16 , taking <p a tt(f) — —£ for r < 
1-5<ths and zero otherwise. By minimizing the grand potential functional we determine the 
equilibrium density profile p(r) and the wall-fluid tension for any thermodynamic state; details 
of the numerical treatment will be given elsewhere. Our DFT approach has the following 
advantageous features: i) the coexisting densities of gas and liquid, p g and pi, can be calculated 
precisely, ii) the DFT obeys the sum-rule © and the Gibbs adsorption theorem and iii) the key 
quantities £ and 7^(00) which enter eqs. <|3I7I8|) can be obtained from independent calculations 
(for the bulk and for the planar interface) using the same functional. 
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Fig. 2 - The difference between contact densities of a square-well fluid in the liquid and in the gas 
phases, at ^t {T), adsorbed at a hard spherical (full line) and hard cylindrical (dotted line) wall of 
radius Ri for fcsT/e = 1. From the linear portion of each curve near the origin, we can extract the 
gas- liquid surface tension 7 9 ;(oo). The next to leading order term is proportional to RJ 2 In R s for the 
sphere and to R~ 2 for the cylinder - see eqs. and 

Fig. 3 - As in fig. H but now for the difference between the surface tensions. For both geometries 
the leading order curvature dependence is proportional to In Ri/Ri - see eq. In both cases the 
coefficient is confirmed to be ai£y g i(oo) / 'UbT ', with £ the correlation length in the bulk gas. 



We focus first on the density profile. In fig. ^ we display results for p{r), at a reduced 
temperature ksT/e — 1 and chemical potential /i^ (T), for the square- well fluid adsorbed on a 
hard sphere with radius R s = 2bQQans- In the same figure we show the profiles for the wall-gas 
interface at /iJ Q (T) and for the planar interface with Sfi = ji — p C o{T) = 27^(00) / \R s (pi — p g ))- 
One can observe that this translation between bulk field for a planar wall and curvature yields 
wall-liquid profiles which lie on top of each other. This is a non-trivial result. Although it is 
evident that within the coarse-grained treatment of eq. the film thickness l l eq is identical 
for planar and curved walls, there is no reason a priori why the complete microscopic density 
profiles should be especially close. 

In fig. |21 we seek to test the validity of eqs. {7| |SJ by examining the difference in the 
contact densities for the two phases as a function of -R" 1 - Our DFT results arc in excellent 
agreement with the prediction of the coarse grained theory for both spheres and cylinders. In 
particular we find the data is fit best with a R~ 2 lni? s contribution for spheres and without 
such a contribution for cylinders. It is straightforward to extract the coefficient of the R~ x 
contribution and we find close agreement, to 1 part in 10 4 , with our independent (planar) 
results for the dimcnsionlcss quantity ai"f g i(oo)aj ls /kBT. 

In fig. we plot the difference between the wall-liquid and wall-gas surface tensions, eval- 
uated at pf a (T), respectively versus Rr l for the square-well fluid adsorbed at spheres and 
cylinders, again for fc^T/e = 1. In accordance with eqs. © and l@J we expect that for the 
sphere the leading order curvature correction to the difference is proportional to Rj 1 \nR s 
with a coefficient given by 2£7 s ;(oo) where £ is the correlation length in the bulk gas. For 
cylinders we employ eqs. @ and JSJ which imply the same type of leading order correction but 
now with a coefficient (£ + bn)j g i (00). Our numerical DFT results are completely consistent 
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with these predictions with the relevant coefficients agreeing with the planar result to better 
than 1%. In our calculations we were able to obtain reliable results for R s up to 30000cths 
and R c up to BOOOcths- For the cylinder our best fit yields a coefficient that has bjj = 
(or a very small fraction of £) and we conclude that there is no evidence for the hiR c /R c 
term conjectured Jlj f° r the cylindrical fluid- fluid interface - see eq. We have confirmed 
that the same level of agreement between the coarse-grained theory and DFT holds at other 
temperatures. 

Since there is nothing special about the square- well model our predictions for hard spheres 
and cylinders should be valid for any fluid with short-ranged fluid-fluid potentials that exhibits 
gas-liquid coexistence and this has important repercussions for the solvation of big solvophobic 
solute particles. The excess chemical potential of a single hard sphere of radius R s in a solvent 
is given generally by fi HS (R s ) — p4nR^/3 + 7^j(i? s )47ri?^, where 7^* is the wall-fluid surface 
tension. Using eq. it follows that the difference between jl HS evaluated in the liquid, where 
drying occurs, and in the gas phase, at fif (T), is 

(r.HS _ ~HS\ 

4^2 = ^(°°)(1 + lni? *) + HOT > ( 9 ) 

for R s — > 00. Whilst it is known from recent studies of hard-sphere solutes that p,f?S 
should contain a gas- liquid surface tension contribution 7 ff ;(oo)47ri?g , the presence of the non- 
analytic correction is striking - especially when we recall that the excess chemical potential 
is the derivative of the excess free energy of the bulk mixture w.r.t. the solute density in the 
limit of vanishing solute. 

The coarse grained treatment that leads to eq. © is valid for wetting by either fluid 
phase. Thus, provided the attractive part of the wall- fluid potential is sufficiently strong and 
has finite range or decays on a length scale that is shorter than the bulk correlation length £, 
now referring to the wetting liquid, adsorption from the saturated gas, at fj,~ a (T), will give rise 
to a wall-gas tension which has a term on £7 9 z( 00 )R^ 1 In Ri. Although the contact theorem 
p)l is modified when the wall-fluid potential is no longer purely hard the contact density in 
the presence of a wetting liquid film will still acquire terms which depend on the gas-liquid 
surface tension, i.e. we expect results equivalent to 10 and ©. 

Both the coarse-grained and DFT approaches presented here are mean-field like in that 
they omit effects of capillary- wave fluctuations in the wetting film For complete wet- 

ting in three dimensions at a planar interface, renormalization group studies based on effec- 
tive interfacial Hamiltonians with the binding potential JQl find that critical exponents are 
not altered from their mean-field values but the amplitude of the equilibrium film thickness 
l eq (oo) is changed from £ to £(1 +u>/2), for u> < 2, when fluctuations are included 0. Here 
u = fcsT/(47T7 g i(oo)£ 2 ) is the usual parameter which measures the strength of capillary-wave 
fluctuations: w — corresponds to mean-field. We conjecture that our present results for the 
leading non-analytic term in the surface tension are modified in a similar fashion, i.e. the 
third term in eq. © simply has £ replaced by £(f + uj/2). As regards the results J7J) and (JSJ) 
for the difference in the contact densities, the leading-order (Laplace pressure) terms will be 
unchanged by fluctuations whereas the next to leading order term in eq. J7J will require the 
same replacement for £. 

For real fluids dispersion forces are always present giving rise to r~ 6 power-law decay of 
the fluid-fluid pair potential. This leads, in turn, to a wetting film thickness which diverges 
at coexistence as l eq (Ri) ~ R^ 3 , for Ri — > 00 ^311201 . The present coarse-grained analysis 
suggests that the non-analyticities in curvature which arise for wetting with such potentials 
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will be power laws rather than terms involving logarithms and we are presently investigating 
these using DFT |2T]. 

We have shown that the isomorphism between bulk field Sfi(pi — p g ) for a planar substrate 
and Laplace pressure at a curved substrate implied by the binding potential leads to 
striking consequences for interfacial properties: first one cannot obtain the surface excess free 
energy (surface tension) of a fluid that wets completely a non-planar substrate by expanding 
only in powers of the curvature (s). Second the true microscopic density profile near the 
curved substrate depends on the surface tension of the gas-liquid interface which, for large 
Ri, can be far from the substrate. A detailed explanation of this curious behavior will be 
given elsewhere. Here we merely state that it is associated with the exponential decay (for 
short-ranged potentials) of the tails of the density profile of the 'free' gas-liquid interface. This 
is not the first time that analysis of complete wetting or drying at a hard wall has caused 
surprises and lead to new insight into the fundamental physics of fluid interfaces |15II22| . 

We have benefited from conversations with J.R. Henderson and M. Thomas. R.E. is 
grateful to S. Dietrich for kind hospitality and to the Humboldt Foundation for support under 
GRO/1072637 during his stay in Stuttgart. 
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